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ABSTRACT 

When a gravitating point mass moves subsonically through a magnetized and isother- 
mal medium, the dynamical structure of the flow is studied far from the mass using a 
perturbation analysis. Analytical solutions for the first-order density and the velocity 
perturbations are presented. Validity of our solutions is restricted to the cases where 
the Alfven velocity in the ambient medium is less than the accretor's velocity. The 
density field is less dense because of the magnetic effects according to the solutions 
and the dynamical friction force becomes lower as the strength of the magnetic field 
increases. 

Key words: magnetohydrodynamics - ISM: general - galaxies: kinematics and dy- 
namics - stars: kinematics 



1 INTRODUCTION 

In many astronomical systems when an object (e.g., a star 
) travels through a collisionless or a collisional system, it 
is decelerated in the direction of its motion because of the 
gravitational interaction of the perturber (i.e, the object) 
and the wakes it excites in the ambient medium. Theoretical 
efforts to study thi s problem have been started after a 
pioneering work by I Chandrasekhaj (| 19431 ). He considered 
a collisionless background and the drag force that the per- 
turber experiences was found analytically. Since then it has 
been applied to many astronomical systems, as in migration 
of protoplanets in protop l anetary accretion dis cs (e.g., 
iMuto. Takeuchi fc IdalliuTll : IVillaver fc Livioll2009l ) , motion 
of supermassive black holes within g as-rich galaxies (e.g., 



Siiacki. Springel fc Haehneit] l201ll : iRuderman fc Spiegell 
197ll). orbita l evolution of embedded binary stars (e.; 



Stahlerl |2010| ). motio n of compact stars around super- 
massive black holes (INaravanl 120001), and the heat i ng of 
the galaxy clusters (IKim. El-Zant fc Kamionkowskil 120051 : 
|E1-Zant. Kim fc Kamionk owski 2004) and even studies at 
the cosmological scales ( Tittlev et al.ll200"H ). 

The standard approach based on the linearized hy- 
drodynamics equations for describing the excited small 
amplitude perturbations and calculating their gravitational 
interaction with the object moving in a straight-line 
traje ct ory (e.g ., Dokuchaevl 1 19641 : iRephaeli fc Salpeterl 
' 1999D . Numerical simulations confirmed 



Il98d : lOstrike 
results of the linear analysis for the drag force (e 
ISanchez-Salcedo fc Brandenburg! 1 19991 : IKim fc Kim] [200' 



The standard approach has be en extende d to the cases 
with cir cular-orbit pertu rber (IKim! [2010) , double per- 
turbers (IKim et alj 2008]) o r accelerated motion in a 
straight-line ( Namounil 20 1 oj ) . In calculating the drag 



force, t he role of the re lativistic effects has also been stud- 
ied by iBaraussel (120071). A mod ified Newtonian approach 
is used by S anchez- S alcedd l|2009l ) to analysis the drag force. 
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Recently, iLee fc Starderl (|201ll ) (hereafter LS) argued 
that the physical size of the perturber is much smaller than 
the accretion radius in the astronomical systems such as 
supermassive black holes within galaxies or giant planets 
inside protoplanetary discs. Under these circumstances, one 
can not neglect the mass accretion and the amplitude of 
the wakes is not small as authors have assumed in their lin- 
ear approach for calculating the drag force. Thus, to a large 
extend it is the mass accretion which provides the linear mo- 
mentum transfer from the ambient medium to the object. 
LS showed analytically that the steady-sate friction force is 
equal to the accretion rate onto the object multiply by its 
velocity. If the accretion rate is prescribed by the Bondi- 
Hoyle rate (|Bondi fc Hovlel["l944l ). the drag force does not 
rise monotonically but instead peaks around Mach number 
0.68 and then begin to decline (LS). It is an interesting re- 
sult that may find application in a number of astronomical 
systems. However, there are many other important physical 
factors that may modify the drag force, and, we think, the 
approach of LS has this possibility to be generalized to in- 
clude those physical ingredients such as a magnetized and/or 
isentr opic ambient medium. Just recently. iKhaienabi fc Dibl 
|2012l ) extended the approach of LS to the case of an isen- 
tropic medium and found a set of analytical solutions for 
the structure of the flow far from object. They proved the 
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proportionality between the friction force and the accretion 
rate also holds in the isentropic case. 

The purpose of this paper is to calculate the dynami- 
cal friction force when the perturber moves with constant 
velocity through a magnetized medium. In most of the as- 
tronomical systems, the ambient medium is actually mag- 
netized, and one may expect the excited wakes are modified 
because of the magnetic effects. In performing this analy- 
sis, we also assume the physical size of the object is much 
smaller than the accretion radius and the approach of LS 
is followed, but including the magnetic fields. Basic equa- 
tions are presented in the next section. Introducing pertur- 
bation expansions in section 3, we will find analytical so- 
lutions for the first-order perturbed variables in section 4. 
However, second-order solutions and the dynamical friction 
force are obtained numerically in section 5. We conclude by 
some physical implications of the results in section 6. 



2 BASIC EQUATIONS 

The basic equation of our problem in the steady state are 
written as 



+ —J x B, 



(1) 

(2) 

(3) 
(4) 



V • ( P U) = 0, 

„(U.V)U = -VP-p^ + l 

V x (U x B) = 0, 
V-B = 0, 

where p, U, P are the density, the velocity and the pressure, 
respectively. Here, R is the radial distance in the spher- 
ical coordinates (R, 9, <p) whose origin is anchored on the 
gravitating mass M (see Figure [TJ). The current density is 
J = V x B. Moreover, we also assume the gas is isother- 
mal, i.e. P = pc\ where c s is the sound speed. We ignore 
the complications from density gradients. Very far from the 
mass which travels in a straight line trajectory with veloc- 
ity V not only the density, but the magnetic fields are both 
assumed to be spatially uniform, i.e. liniH^oo p = po and 
lmiij-joo B = Boe z where e z is the unit vector along the z 
direction (Figure [T]). 

Introducing a set of non-dimensional variables as r = 
R/Rs, u = U/cs, g = p/po and b = B/£>o, we can write 
our main equations in the non-dimensional forms. Here, the 
sonic radius is defined as R B = GM/c 2 . Therefore, we have 



V • (gu) = 0, 
g(u- V)u = -X7g 

V x (u x b) = 0, 
V-b = 0, 



4+C 2 (Vxb)xb, 



(5) 

(6) 

(7) 
(8) 



where the non-dimensional parameter £ is defined as the ra- 
tio of the Alfven velocity Va = Bo/V^rpo and the sound 
speed, i.e. £ = Va/c b . Obviously, in equations ©-([8]) the 
radial derivatives are calculated with the respect to the non- 
dimensional radial distance r. We also define j3 = V/cs. 
Thus, properties of the flow at the distances far from the 
gravitating object M are u = /3e z , g = 1 and b = e z when 
r » 1. 



r= R/R 



V 




B, 



M 



Figure 1. This figure shows the assumed magnetic field and the 
flow velocity. The object is fixed at the origin and the magnetic 
field and the velocity of the gas are both uniform at the large 
distances from the object. 



Note that in our subsequent calculations, it is assumed 
that the system is axisymmetric (d/dtp = 0). If we introduce 
the non-dimensional stream function ip(r,9) as 



gu = V x 



(9) 



v r sin 6 

then the non-dimensional continuity equation ([5]) is auto- 
matically satisfied. Thus, the components of the velocity can 
be written in terms of the stream function as 
1 dtp 



Ug = 



gr 2 sin 6 89 



dr 



(10) 



(11) 



gr sm f 

We can also define the non-dimensional positive-definite 
magnetic flux function $ as 

<h 



b = V x (- 



(12) 



So, equation ((HJ is satisfied and the components of the mag- 
netic field become 



6 r = 



r 2 sin 9 89 

-l a$ 

■ sin 9 dr 



(13) 



(14) 



Having the above equations for the components of the 
velocity in terms of ip and the magnetic field in terms of "3>, 
we are left only with the momentum equation ((6| and the 
induction equation 0. Thus, we have 



dr 



Ug dll r 



d . . db T ~ 
d-r {rb0) '^9. 

I dug Ug dug 



U^Ug 



d_ 

dr 
u r bg 



rbg) - 
- Ugbr 



r d9 


db r - 




d9 _ 




= 0. 





dg g c 2bg 
or r z r 



19p 2&r 

\d9 + t r X 



(15) 

(16) 
(17) 
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Figure 2. Isodensity curves corresponding to the first-order den- 
sity perturbation for = 0.5 and £ = 0.4 (solid) and £ = 
(dashed). Each contour is labeled by its level. The gray inner cir- 
cle denotes the sonic radius. Obviously, the solutions are more 
accurate at the distances far from the inner sonic circle. 



Since we have neglected the dissipative terms in the induc- 
tion equation and the system is steady state, we could say 
based on the equations ([5]) and ([8]) that the magnetic field 
and the velocity vectors are parallel. The simplified induc- 
tion equation (jXTJ) also shows this behavior. 



3 PERTURBATION EXPANSIONS 

We expand the density and the stream function similar to 
LS as 

g(r, 9) = 1 + MtflK 1 + g-2{9)r- 2 + g- 3 (6)r- 3 + • ■ • , (18) 

i>{r,9) = f 2 (9)r 2 + h(6)r + f (e) + /-iWr" 1 + ■ • ■ , (19) 

and in our magnetized model, the magnetic flux function is 
expanded as 



$(r, 6») = h 2 (9)r 2 + h^{6)r + h {9) + /i_i(0)r 



+ 



(20) 



Since the velocity tends to be constant far from the 
object, one can simply obtain f 2 {9) = (l/2)/3 sin 2 9. Also, 
magnetic field lines are uniform far from the object and then, 
h 2 (9) = (V 2 ) sin2 e - We can set </>0:7i") = 0, and so f,(n) = 
for i = 1,0, —1, —2 etc., and regularity of the components 
of the velocity on the axis defined by 9 = and 9 = n 
implies fl(ir) = /i(0) = 0, for i = 1,0,-1,-2, etc., and 
/i(0) = for i = 1, —1, —2, etc. (see LS for details of their 
arguments). Boundary conditions on the geometry of the 
magnetic field lines also imply hi (it) = h'^n) — h[(Q) — 0, 
for i = 1,0, —1, —2, etc., and hi(0) = for i = 1, —1, —2, etc. 
Substituting the above expansions into the main equations, 
a set of ordinary differential equations are obtained, which 
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Figure 3. Profiles of fo(0) and g~2(8) for different values of £ 
when the initial value of g_2(7r) is -0.88. Each curve is labeled by 
the corresponding value of £. 



are solvable analytically and/or numerically subject to the 
mentioned boundary conditions. 

Also, the mass accretion rate is written as 



M = -2tt / pU n R sin 9d9. 
/o 



(21) 



If we nondimensionalize the accretion rate M by 2irpoc B R 2 
and using expansion 1)191) . the nondimensonal accretion rate 
M becomes (see LS for the details) 



M = /o(0). 



4 FIRST ORDER EQUATIONS 



(22) 



We can now substitute the expansions (|18|) . (|19[) and (|20l) 
into the momentum equations (|15|) and ()16[) and the induc- 
tion equation (|17[) . Then, we match the coefficients of each 
power of r. Like the non- magnetized case (LS), we find the 
highest power of r is r , but all the corresponding coeffi- 
cients on both sides of the equations vanish identically. So, 
we consider the next power of r (i.e. r~ 2 ) and the first-order 
equations are obtained. From the radial component of the 
momentum equation, we have 



-Pf" ~ Pfi + P sin 9 cos 
1 + € 2 (^i - cot 9h[) = 0, 



+ (P 2 cos 2 9- 1)5-1+ 



(23) 
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g ,00 = 0.12 

P = 0.5 _ 




0.0 0.2 0.4 0.6 0.8 1.0 

0/ti 

Figure 4. The same as Figure[3] but for g_2(7r) = 0.12. 

and the 9— component of the momentum equation gives 
(1 — p 2 sin 2 9)g'_ 1 — p 2 sin 8 cos 9g~i + 

i 2 cot 9(h" - cot 9h[ ) = 0, (24) 

and finally the induction equation (|17[1 yields 

cos 9fi - P cos 9hi + P sin 9h[ - sin 9f{ = (25) 

For the non-magnetized gas, the above first-order equa- 
tions reduce to the equations of LS if we set £ = 0. We can 
solve these equations for fi, g~i and hi, and then proper- 
ties of the flow to the first order are determined. One can 
integrate equation (|25[) simply by re-arranging its terms as 

cos0(/i - Phi) - sin6>4(/i - ph x ) = 0. (26) 
a9 

Introducing W = fi — Phi, this equation becomes dW/d9 — 
W cot 9 and its solution is W = Wo sin 9 where Wo is an 
arbitrary constant to be determined from the boundary 
conditions. Thus, equation /i — ph\ = Wo sin 9 is valid 
for the whole range of 8 including 9 = n. So, we have 
/((•7r) — ph'i(n) = Wjcos7t = — Wo. Considering our im- 
posed boundary conditions /{(tt) = ^(tv) — 0, we obtain 
Wo = and 



(27) 



fx = Phi. 

We can now substitute equation (I27|l into equations (|23p 
and (|24p and the following equations for /i and g~i arc 
obtained 




0.0 0.2 0.4 0.6 0.8 1.0 

Q/n 

Figure 5. The same as Figure[3] but for g_2(7r) = 1.12. 



~ j 0ft - Pfi + P 2 sine cos 9gU+ 



(P cos 



l) ff _i+l = 0, 



e 2 

(1 - P 2 sin 2 9)g'_ 1 - p 2 sin 9 cos 9g- X + 5- cot Of ['- 



^cot 2 0f[ 



0. 



(28) 



(29) 



Thus, equations ()28|) and (|29|) are our main equa- 
tions to be solved. Interestingly, these equations are inte- 
grable, though the mathematical manipulation is cumber- 
some. Here, we summarize the solutions as 



m(«) 



i - e/P 2 



VT - p 2 sin 2 



(30) 



(31) 



where T = 1 + £ 2 - £ 2 / P' 2 ■ For the non-magnetized case 
(i.e. £ = 0), we have T = 1 and the above solutions reduce 
to equations (25) and (26) of LS. Note that the above first 
order solutions are physically acceptable if T — p 2 sin 2 9 > 0. 
If for 9 = 7r/2 (where the second term has a maximum value) 
the left hand side becomes positive, then for other values of 
9 the left-hand side would be positive (as we want). Having 
r = 1 + £ 2 - £, 2 /P 2 , we can write 1 + £ 2 - ( 2 j p 2 - p 2 > 
which implies /3 > £. Thus, our solutions are applicable to 
the case where the Alfven speed is less than the sound speed. 
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Figure 6. Density contours for the non-magnetized (top plot) 
and the magnetized flow (bottom plot) with /3 = 0.5, § = 0.4 and 
g_ 2 (7r) = -0.88. 



Figure [2] shows isodensity contours for the first-order 
density perturbation, i.e. equation (131[) . Solid contours are 
for a magnetized ambient medium when £ = 0.4. Nonmag- 
netized contours are also shown by the dashed curves. In 
both cases, the Mach number is fixed /3 = 0.5. A circle with 
unit radius at the center denotes the sonic radius. Obvi- 
ously, our solutions are valid beyond the sonic radius and 
all the contours are drawn for the radii larger than 2. Each 
contour is labeled by the corresponding level. For the non- 
magnetized solutions, the density varies from 1.09 to 1.58 at 
the inner parts and the magnetic effects change this interval 
from 1.09 to 1.48. In other words, the density field in the 
presence of the magnetic field is less dense in comparison to 
the non-magnetic solutions. When magnetic fields are not 
considered, the perturbations in the ambient medium grow 
easier and so the wakes have larger amplitudes comparing 
to the magnetic case. It is because of the magnetic pressure 
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Figure 7. Streamlines for the non-magnetized (top plot) and the 
magnetized flow (bottom plot) with /3 = 0.5, £ = 0.4 and g—2(ir) = 
-0.88. 

which acts as an extra pressure against further growth of 
the perturbations. 



5 SECOND-ORDER EQUATIONS 

Now, we can equate coefficients of the next power of r in the 
components of the momentum equation and the induction 
equation. From the induction equation, we obtain 

fo = /3h . (32) 

Considering the above relation between /o and ho, we 
obtain the following equation from the radial component of 
the momentum equation 

(j ~ P)fS - W + j) cot ef + f sin cos 8g'_ 2 + 

(2f3 2 cos 2 9 - 2)ff-2 = Ai + Ai + Az + Ai, (33) 
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where all the terms on the righthand side depend on the first 
order variables as 



f? httcosO , (/0 2 , htt 



+ 



+ ■ 



(34) 



A2 = Pfig-i - 2/3/( <?_, cot 9 - PfigU cot 9- 

Pfig'-i+Pfi'g-i, (35) 

_4 3 = 2(< ? _ 1 ) 2 -3< ? _ 1 . (36) 
p p sin 6 



fl'f 



P sin 2 9' 

Also, the 6>— component of equation of motion gives 

(£) cot 0/o - (/3 + ^ cot 2 0)/g + (1 - p 2 sin 2 9)gL 2 - 

2p 2 sm6cos6g- 2 = Bi + B 2 + B 3 + B 4 , 
where 

_ 2 cotfl /x/{ 
1 71 sin 2 6 sm 2 0' 
B 2 =Phg-i cot 9 + Pfig-! + 2Pf 1 g'_ 1 , 
B3 = -2fl-iffli, 

B < = ^ - 4fi - Wi'g-i cot 0+ 

p^ sin 9 sin # 
2pf[g. 1 cot 2 9). 



(37) 



(38) 

(39) 

(40) 
(41) 

(42) 



Both righthand sides of equations (|33[) and ()38[) are de- 
termined analytically because we have already found /1 and 
g~i, i.e. equations (|30[) and (1311) . Although equation (|38l) is 
integrable as we will show, it is very unlikely to solve this 
equation and equation (|33[) for fo(9) and g~ 2 (9) analytically. 
So, one should solve these equations numerically subject to 
the appropriate boundary conditions. Note that if we set 
£ = 0, these equations reduce to equations (27) and (31) 
of LS. Introducing D = F — p 2 sin 2 9, we can re-write the 
second-order equation (|33[) as 



(j ~ P)fo ~ (P + j) cot 0f o + P 2 sin 9 cos 9gL 2 + 
(2P 2 cos 2 9 - 2)g_ 2 ' 



-r(r-o)D 2 (/3 2 - l) 2 

(7li£> 5/2 + Tl 2 V 2 + TZ 3 V 3/2 + TliV + TZ 5 ), (43) 
where 

TZj = VT(1- P 2 )(P 2 +2-3T), (44) 

TZ 2 = 2F 3 - 4/3 2 r 2 + (5P 2 - 3)r + p 2 - /3 4 , (45) 

ft 3 = 3r 3/2 (/3 2 - l)(/3 2 - T), (46) 

Hi = r(r - p 2 )[-4r 2 + 2(1 + 2p 2 )r + 1 - 3p 2 ], (47) 

^ 5 = 2r 2 (r- i)(p 2 -r) 2 . (48) 

Also, equation (|38|) is written as 



(L) cot 0fH - (p + L cot 2 9)f + (1-P 2 sin 2 9)gL 2 - 



2/3 sin # cos ( 



P 2 sin 61 cos 9 



-r(r -V) 2 V 2 {P 2 - l) 2 
(ft 6 £> 5/2 + II7V 2 + TZsV 3/2 + TZ 9 V + 1Z W V 1/2 + Ku), (49) 
where 

ft 6 = \/r(r-i)(/3 2 -i), (50) 
^ 7 = 2r 3 - 4p 2 r 2 + {5p 2 - 3)r + p 2 - p"\ (51) 

^ 8 =r 3/2 (4/3 2 -l-3r)(/3 2 -l), (52) 



tz 9 = r(r - /? 2 )[-4r 2 + (5p 2 - i)r + 2 - 2p 2 



(53) 



H10 = 2r 5/2 (r-/? 2 )(/3 2 -i), (54) 
= r 2 (^ 2 - r)[-2r 2 + (3/? 2 - i)r + 1 - p 2 ]. (55) 



The lefthand side of the second-order equation (|49p is 
perfect derivative. By integrating both side of this equation 
from 9 — n to 0, we find 



f 2 f 2 
5Lcot0/5 + (5L 



p)f +Vg- 2 



Since /o(tt) = ,fo(0) = 0, the above equation gives 



/o(0) = 



P 2 -t 



■fa- a (0)- S _ 2 (ir)]. 



(56) 



(57) 



If we set £ = 0, the above equation reduces to equation (46) 
of LS. Here, like to the non-magnetized case, the accretion 
rate is proportional to the difference, upstream and down- 
stream, of the second-order density perturbation. However, 
the proportionality constant not only depends on the veloc- 
ity of the accretor, but also it depends on the magnetic field 
strength. 

We can solve the equations for the same three values 
of g- 2 (n) as were used by LS, just in order to make easier 
comparison to the non-magnetized case. Note that our main 
equations Q43p and (|49|) are singular at both endpoints 9 = 
and 9 — n. Figures [3] [4] and [5] show behaviors of fo(9) and 
g- 2 {9) for three fiducial values of g- 2 {ir) = -0.88, 0.12, 1.12. 
Just to make an easier comparison to the non-magnetized 
profiles of LS these boundary values for g- 2 (n) have been 
selected. In each figure, the starting value of g- 2 (n) is fixed, 
and we play with the ratio £ as a free input parameter. 
Profiles of the second-order coefficients strongly depend on 
the strength of the magnetic field as shown in these figures. 
The non-magnetic curves passes through the same point at 
9 — n/2. But magnetic effects change this trend, in partic- 
ular, when £ becomes closer to p\ However, the coefficient 
g- 2 {9) is flat at 9 = and n which means the second-order 
density perturbation is uniform at both the downstream 
and upstream of the flow irrespective of the existence of 
the magnetic fields. It is understandable based on our im- 
posed boundary conditions. Density contours and the veloc- 
ity streamlines are also shown in Figures [6] and [7] for the 
magnetized and the non-magnetized case corresponding to 
the parameters P = 0.5, £ = 0.4 and g- 2 (ir) = —0.88. Here, 
the curves are shown for r 2. But the inner contours in 
these figures continue to the region with r < 2 which is not 
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Figure 8. The nondimensional dynamical friction force versus 
the Mach number /3 for different values of £. Each curve is limited 
by the inequality £ < and the numbers are the corresponding 
value of £. 



shown and so, are connected vertically artificially. The rea- 
son we exclude an inner part goes to back to the validity 
of our expansions for the radial distances greater than the 
accretion radius. 

Like LS, there is no restriction on <?_2(7r) in the mag- 
netized case as long as we do not extend the solutions to 
the inner part of flow or else adopting a model for the ac- 
cretion and fixing the accretion rate. LS followed the sec- 
ond approach which is much easier and avoids complexi- 
ties related to the conditions on the sonic s urface. Their 
accretion model is t he classical Bondi-Hoyle (|Bondil I"l952t 
iBondi fc Hoviell 19441 ) model for which there is an analytical 
relation for accretion rate as, 



M 



2(A 2 + /3 



2x1/2 



(1 + /3 2 ) 2 



(58) 



Having the above the relation for the accretion rate, one 
boundary condition is given at 9 — as /o(0) = M . The rest 
of boundary conditions are at 9 = n as /o(tt) = fa{it) = 0. 
Thus, we can solve the second-order equations as a two- 
points boundary value problem. As far as we know, for the 
magnetized accretion there is not a closed analytical for- 
mula for the accretion rate. However, it was shown that 
the magnetic field effects suppress the accretion rate (e.g. , 
Igumenshchev 20061; IShadmehril |2004|; IShcherbakovl [2008; 



Igumenshchev k. Naravanll2002l ; IPerna et al.ll2003h . For ex- 
ample, when the accretion occurs as a radiatively inefficient 
flow, the magnetic field is amplified as the gas flows in, and 
so, the accretion rate is reduced i n comparison to the Bondi 
rate a s studied numerically by llgumenshchev fc NaravarJ 
(2002). But there are many uncertainties regarding to the 
magnetized accretion, in particular about the true physical 
mechanism of the accretion or when the accretor is in trav- 
eling. In order to proceed analytically, we introduce a toy 
accretion rate as 



JWmag - ,j + £ 2 + , (OV) 

which gives a reduced rate because of the magnetic field. 
Here, the rate decreases with £ and it gives Bondi-Hoyle ac- 
cretion rate if we set £ = 0. It is as if we replace the thermal 
pressure by the sum of thermal and the magnetic pressures, 
i.e. c 2 — > c 2 + v\. In fact, the extra pressure provided by 
the magnetic force lead to the reduction in the accretion 
rate. We emphasize our prescribed accretion rate relation 
(|59[1 is not an accurately physically confirmed relation, but 
it shows a reduction to the accretion rate with increasing 
the strength of the magnetic field at least qualitatively. 

Having equation (I59|l for the accretion rate, we did a 
detailed parameter study by solving the main equations (|43|l 
and (|49[) numerically. We also obtained the non-magnetic 
results of LS by setting £ = 0. In particular, LS proved 
that the nondimensional dynamical friction force for a non- 
magnetized flow is 



F 



[(1 - /3 2 ) sin 9 cos 9g- 2 + /?(! + cos 2 9)fg] d9, (60) 



and then they showed analytically that the above integral 
is equal to f3M. They obtained the above relation by inte- 
grating the stress tensor over sphere surrounding the object. 
Since we are considering magnetic effects, it is known that 
the contribution of the magnetic field to the stress tensor is 
^B 2 Sij — BiBj, i.e. the tress tensor is 

T = pUU+(P+is 2 )I-BB. (61) 

Thus, in our non-dimensional notation, the magnetic part 
of the drag force becomes 



F lr 



/ /) /■" Ci IS (I Sll] I 



1 + t 2 J bzbrr 2 sin6d8. (62) 



Also, we have 



, , 2 . „ cot6» /9$\ 2 1 

bArSm9 = —{d0) + r(W{^) 



and 



r 2 " {'89 ) + V"cV ) 



(63) 



(64) 



Now, we can substitute corresponding expansions into 
the above equations to calculate the drag force. However, 
the resulting integral is not integrable and one should cal- 
culate it numerically. Our numerical solutions show that the 
dynamical friction force in the magnetized case is approx- 
imately equal to (/3 — £ 2 //3)/o(0), where the accretion rate 
is replaced by the magnetized accretion rate, i.e. equation 
((59)) . Thus, we numerically confirm F ~ (/3 - (, 2 //?)M mag . 
Figure [S] shows dynamical friction force F versus parameter 
/3 for different values of £. Note that our solutions are re- 
stricted to the cases with £ < /3. Each curve is labeled by its 
parameter £ and it is truncated due to the inequality £ < j3. 



6 CONCLUSIONS 

When a gravitating object move subsonically through a 
magnetized medium, we studied the structure of the flow 
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using a perturbation analysis. Our study is a direct gener- 
alization of LS to the magnetized case. We found analytical 
relations for the first-order density and velocity perturba- 
tions. According to these solutions, the structure of the flow 
in comparison to the non-magnetized medium is less dense. 
However, in order to determine the dynamical friction force 
the second-order variables are needed. The dynamical fric- 
tion force was found to be F ~ /3M mag which implies a 
reduced force because of considering the magnetic effects. Is 
this relation valid for a physically confirmed magnetized ac- 
cretion rate? This is an open question to be studied further. 
Future models of the magnetized accretion should give us 
an analytical or numerical relation for the accretion rate as 
a function of /3 and £ , i.e. M mag (/3, £). Then, our main equa- 
tions (|43p and (|49[l could be solved easily to calculate the 
dynamical friction force. However, we think the magnetic 
effects will still reduce the dynamical friction force. 

In the cluster of the galaxies, the motion of 
the each member of group is affected by the dy- 
namical friction force and it has been suggested 
that work done by the force may heat up the in- 
traga la ctic medium llEl-Zant. Kim fc Kamionkowskil 



120041 : iKim. El-Zant fc Kamionkowskil 120051 ). This ex- 
tra source of the heating may resolve the cooling flow 
problem in the cluster of the galaxies (for a detail ed 



analysis, see IKim. El-Zant fc Kamionkowskil 120051) 1 



IKim. El-Zant fc Kamionkowskil (|2005l ) used the classical 



formula for calculating the drag force which is based 
on assumption that the object's radius is greater than 
the accretion radius. But as LS argued in such systems, 
accretion radius is comparable to the radius of the accretor. 
Thus, it would be interesting to calculate the heating due 
to the friction force using the formula that is presented 
by LS, at least to the regions of the cluster where the 
galaxies are moving nearly subsonically like at the outer 
parts. But intracluster medium is magnetized, and so, our 
analysis imply a lower heating rate when magnetic fields 
are taken into account. We finally note that several issues 
of potential importance were not investigated in this paper. 
For example, we assumed that the ambient medium is 
isothermal. One can relax this assumption and re-do a 
similar analysis, but for a magnetized isentropic ambient 
medium. 
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